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t^ - '. Abstract 

The magnetic helicity has paramount significance in nonlinear saturation 
of galactic dynamo. We argue that the magnetic helicity conservation is 



violated at the lepton stage in the evolution of early Universe. As a result, 
a cosmological magnetic field which can be a seed for the galactic dynamo 
obtains from the beginning a substantial magnetic helicity which has to be 
O . taken into account in the magnetic helicity balance at the later stage of 

galactic dynamo. 

,§ : 1 Introduction 

CO 

a 

Magnetic fields of galaxies are believed to be generated by a galactic dynamo 
based on the joint action of the so-called a-effect and differential rotation. 
The a-effect is connected with a violation of mirror symmetry in MHD- 
turbulence and therefore caused by rotation. For a weak galactic magnetic 
field, the mirror asymmetry is associated with helicity of the velocity field 
and is proportional to the linkage of vortex lines. 

The magnetic helicity is an inviscid integral of motion and its conserva- 
tion strongly constrains the nonlinear evolution of the galactic magnetic field. 
The helicity density B • A of a galactic large-scale magnetic field is enhanced 
by galactic dynamo (here B is a large-scale magnetic field, A is its vector- 
potential). Because of magnetic helicity conservation, this income must be 
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compensated by magnetic helicity of a small-scale magnetic field. Note that 
the magnetic helicity density is bounded from above by bfl, where bi is the 
magnetic field strength at the scale /. Hence the supply of a small-scale mag- 
netic helicity occurs to be insufficient for the compensation required. This 
fact strongly constrains the galactic dynamo action (see [1, 2] and references 
therein) . 

On the other hand, a weak almost homogeneous cosmological magnetic 
field B c can introduce a new element in this scheme. If we adopt a thick- 
ness of the gaseous galactic disc l ga \ ~ 1 kpc as a typical scale for galactic 
dynamo action and B c = 10~ 9 G, then the magnetic helicity supplied via 
the cosmological magnetic field (h c ~ B 2 Ih ~ 10 10 G 2 cm, where Ih is the 
horizon size) can be of the same order as that one for the galactic magnetic 
field (/i gal ~ B G l gal ~ 3 x 10 9 G 2 cm, B G « 1CT 6 G). 

The cosmological magnetic field B c if it exists must be substantially 
weaker than the galactic magnetic field (see [3] for review). According to 
the analysis of rotation measures of remote radio sources, B c < io -10 - -11 G. 
This estimate however is based on the assumption that a substantial part 
of charged particles in the Universe is in form of thermal electrons in the 
intergalactic medium. A more robust estimate, B c < 1CT 9 G is based on the 
isotropy constrains. 

We estimated above that if the cosmological magnetic field is of the order 
B c = 10~ 9 G its magnetic helicity density can be comparable with the mag- 
netic helicity density of the galactic magnetic field. Of course, the estimate 
B c = 1CT 11 gives a much lower value for the magnetic helicity density. How- 
ever it is more than natural to expect a magnetic helicity concentration in 
the processes of galactic formation. The question is whether a mechanism of 
magnetic helicity production can be suggested for physical processes in the 
early Universe. 

Here we suggest a mechanism for magnetic helicity generation by a collec- 
tive neutrino-plasma interactions in the early Universe after the electroweak 
phase transition. 
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2 Weak force contribution to dynamo action 
in early Universe 



Let us consider the electron-positron plasma as a two-component medium, 
for which the strong correlation between opposite charges due to Coulomb 
forces gives V± ~ V. Here V is the common fluid velocity of electroneutral 
conducting gas while its positively and negatively charged components have 
different velocities, V± = V±5V, where a small difference <5V <C V gives the 
separation of charges at small scales and enters the electromagnetic current 
j(em) = 2 | e | n e 7 e 5V obeying in MHD the Maxwell equation, 47rj( cm ) = 
VxB(fi = c=l). 

The electric field E derived from Euler equations for plasma components 
takes the form [4] which includes the contribution of weak interactions E weak 
taken in the collisionless Vlasov approximation. We do not consider other 
known terms which describe weak interaction collisions [5] , Biermann battery 
effects, etc., and which do not play substantially in favor of helicity genera- 
tion. In turn, we keep in E wca k only the axial vector term which violates the 
parity: 



is the axial weak coupling, upper (lower) sign is for electron (muon or tau) 
neutrinos, 5n Ua = n Ua — n-pjs the neutrino density asymmetry, 5j Ua = j„ a —jv a 
is the neutrino current asymmetry; b is the unit vector along the mean 
magnetic field; no is the lepton number density at the main Landau level 
given by the equilibrium Fermi distribution in a hot plasma, n ~' ) ~ n + ^ = 
(| e | -B/27r 2 )Thi2, where at equilibrium the temperature obeys T + = T_ = 
T ^> m e . In the non-relativistic limit, v <C 1, the relativistic polarization 
terms in Eq. (I) tend to the lepton densities at the main Landau level, 



Accounting for the first line in Eq. (I) we obtain the axial vector term 
E^ ak = -«B, where the helicity coefficient a is the scalar in the standard 
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model (SM) with neutrinos instead of the pseudoscalar (v ■ (V x v))) in 
standard MHD: 



G 

a 



2V2\e\B^ 



n P I dt 



In 2 / 1Q- 5 T \ fSn u 



4^ 



(2) 



Here we substituted n u /n e = 0.5, and assumed a scale of neutrino fluid 
inhomogeneity t ~ A{^ id , that is small compared with a large A-scale of the 
mean magnetic field. 

Let us stress that instead of the difference of electron and positron con- 
tributions in axial vector terms entering the pair motion equation [4] and 
given by the polarized density asymmetries ~ (uq ^ — n^) we obtained here 
the sum of them ~ (n^ + n^) that can lead to a significant effect in a hot 
plasma. 

The admixture of the pseudovector «B to the pure vector E, e.g. for 
the constitutive relations D = eE + /3B, H = 7E + /i _1 B due to the same 
neutrino-plasma weak interaction described by the constants f3, 7, has been 
already discussed in literature (see [6], Eqs. (3.5), (3.6)). In a forthcoming 
paper we show that such unusual coefficients appear in chiral media and are 
simply connected with a given by Eq. (2), (e — l)a — (5 + 7, where e is the 
dielectric permittivity of plasma. 

Thus, using Eq. (1) from the Maxwell equation <9 t B = —V x E one 
obtains the Faraday equation generalized in SM with neutrinos and antineu- 
trinos: 

— = V x aB + r]V 2 B , (3) 

where we omitted the weak vector contribution V x E^ k ~ V x <5ji, (x, t) 
suggested by Brizard et al (2000) since we neglect any neutrino flux vorticity 
in the hot plasma of early universe. Because the early universe is almost 
perfectly isotropisc and homogeneous, we ignore here any contribution from 
large-scale motions as well. In the relativistic plasma the diffusion coefficient 
rj takes the form rj = (4tt x 137 T)" 1 . 
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The first term in r.h.s. of Eq. (3), V x «B, is associated with the parity 
violation in weak interactions in the early universe plasma. 

We stress that the Eq. (3) is the usual equation for mean magnetic field 
evolution with a-effect based on particle effects rather on the averaging of 
turbulent pulsations. It is well-known (see e.g. [7]) that Eq. (3) describes 
a self-excitation of a magnetic field with the spatial scale A m rj/a and the 
growth rate a 2 /4r]. Authors [8] estimated these values for the early universe 
to get 



Here l H (T) = (2H)- 1 is the horizon size and H=4AQ x l(r 22 (T/MeV) 2 MeV 
is the Hubble parameter; the variable x = T/2 ■ 10 4 MeV corresponds to 
the maximum temperature T ~ 20 GeV <C T EW ~ 100 GeV for which the 
point-like Fermi approximation for weak interactions we rely on is still valid. 
Finally l u {T) is the neutrino free path and ^ Ve = fi Ue /T <C 1 is the small 
dimensionless electron neutrino chemical potential normalized in (5) on the 
maximum value £„ e < 0.07 allowed by the Big Bang Nucleosynthesis (BBN) 
bound on light element abundance [9]. 

Thus, while in the temperature region T E w > T » T = 10 2 MeV 
there are many small random magnetic field domains, a weak mean magnetic 
field turns out to be developed into the uniform global magnetic field at 
temperatures below T (see Eq. (4). The global magnetic field can be small 
enough to preserve the observed isotropy of the cosmological model [10] while 
being strong enough to be interesting as a seed for galactic magnetic fields. 
This scenario was extensively discussed by experts in galactic magnetism 
[11], however until now no viable origin for the global magnetic field has 
been suggested. We believe that the dynamo based on the a-effect induced 
by particle physics solves this fundamental problem and opens a new and 
important option in galactic magnetism. 
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3 Magnetic helicity generation by collective 
neutrino-plasma interactions 

Let us consider how the collisionless neutrino interaction with charged leptons 
can produce the primordial magnetic helicity H — J V (A ■ H)d 3 x, where v is 
the volume that encloses the magnetic field lines. 
For that we should substitute into the derivative, 



the electric field E given by Eq. (1). Neglecting any rotation of primordial 
plasma given by the first dynamo term, or retaining the resistive term and 
the weak interaction term E^ k given by Eq. (1) that is the main one in the 
absence of any vorticities, one finds from (6) 



Note that the second term in the r.h.s. violates parity: it is a pure scalar 
while other terms are pseudoscalars as it should be for the helicity H in 
standard MHD. Nevertheless, all terms in the generalized helicity evolution 
equation (7) obey CP-invariance as it should be for the electroweak interac- 
tions in SM since the new coefficient a (2) is CP-odd, (CP) a (CP) -1 = —a, 
as well as V x .... This is due to the changes n _ < — > n 0+ and 5n Ua — > — Sn Va 
in (2), provided by the well-known properties: particle helicities are P-odd 
and particles become antiparticles under the charge conjugation operation C, 
in particular, active left-handed neutrinos convert to the active right-handed 
neutrinos under CP-operation, v a — > V a . Obviously, the product (E ■ B) 
entering the helicity evolution (6) is CP-odd too because both electric and 
magnetic fields are C-odd and have opposite P-parities. 

First term in the r.h.s. of Eq. (7) gives conventional ohmic losses for 
magnetic helicity and usually is neglected in helicity balance. 

4 Seed magnetic helicity in cosmology 

Neglecting in evolution equation (7) the first diffusion term, we can calculate 
the magnetic helicity H(t) using Eq. (5) to yield 




(6) 
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dt 



H{t) = 2B 2 max Jd 3 rJdt'a(t>) >J + tf(t max ), (8) 

^ ^max 

where H(t max ) is the initial helicity value at the moment t max if it ex- 
ists, and we present the helicity density entering the integrand as h(x) = 
2.4 x 10 3 (B max ) 2 m~ 1 J(x). Here we use the dimensionless variable x = 
T/2 ■ 10 4 MeV. The maximum value x = 1 corresponds to the maximum 
temperature T max ~ 20 GeV <C T EW ~ 100 GeV (see motivation above), and 
we assumed that the neutrino gas inhomogeneity scale is of the order of the 
neutrino free path l„(T), X^j id ~ l v - The WKB value of the mean magnetic 
field amplitude B max obeys B max < T^ ax /e = 4(T max /MeV) 2 B c , where the 
Schwinger field B c = m 2 e /e = 4.41 x 10 13 G. Hence we may introduce a small 
WKB parameter k <C 1 for the mean field -B max = nT^Je <C T^ w /e, or sub- 
stituting T max = 20 GeV one obtains B max = kx 7x 10 22 G < B EW ~ 10 24 G, 
where, let us say, k ~ 0.01. 

It is worth noting that such WKB value of -B max which is scaled be- 
ing frozen-in as B(T) = _E> max (T /T max ) 2 does obey the BBN limit B < 

10 11 G at the temperature T BBN ^ 0.1 MeV, i.e. B(T BBN ) = k x (7/4) x 

10 12 G. Note also that the sign of the first term in (8) is not well deter- 
mined since it depends on the combined neutrino density asymmetry [8], 
Sn„/n v = Y.a c w ) J n vJ n va ~ &v + _ &d» where the values of the di- 
mensionless neutrino chemical potentials £ Va = fi Ua /T are given by the BBN 
limit citeDolgovl : -0.01 < f„ e < 0.07, or by the CMBR/LSS bound [12]: 
—0.01 < £„ e < 0.22, I £ VfitT |< 2.6. One can use, e.g., the conservation 
of the lepton number L e — L M that implies £„ e = —£ v , however, this does 
not guarantee the definite sign of the combined neutrino density asymmetry 
5n u /n u . The definite sign of the magnetic helicity (left-handed, H < 0), 
arising during electroweak baryogenesis [13] is another case connected with 
the CP-violation. 

Let us emphasize that cosmological helicity production via the collective 
neutrino interaction with hot plasma ceases if the neutrino chemical potential 
fi Ue (hence the neutrino density asymmetry 5n u /n u ) vanishes, £„ e = H V JT — > 
0, 5n u /n v — > 0. Solely the inequality £ Ve < 0.07 is known from the BBN 
bound on light elements abundance at T < O(MeV) (Dolgov et al, 2002), 
thereby substituting for a rough estimate £„ e = 0.07 we estimate the integral 
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for h(x) as J(x) ~ 10. 

Thus, collecting numbers for -B max , J(x) and using the electron Compton 
length m~ l = 3.86 x l(T n cm, one finds the huge value of cosmological 
helicity density that could seed galactic magnetic helicity 

h(x) ~ 4.5 x 10 3 V J(x) G 2 cm ~ 4.5 x 10 3 V G 2 cm. (9) 

5 Discussion 

Traditional galactic dynamo considered galactic magnetic field produced from 
a very weak seed field. This implies that the magnetic helicity of the seed 
field is weak. We argue that the applicability of this viewpoint is limited. 
The seed field for galactic dynamo can be a field of substantial strength and 
substantial helicity. The first part of this statement is already quite well- 
accepted in modern galactic dynamo (see e.g. [3]) while the second one is 
new. In this letter we suggest a physical mechanism for the magnetic helicity 
production for the seed field of galactic dynamo. As far as we know, such 
mechanisms was not considered previously. Note that [14] pointed out in a 
general form the importance of the electroweak phase transition for magnetic 
helicity generation. 

We stress that the epoch just after the electro-weak phase transition and 
that one of galaxy formation are quite remote in respect to their physical 
properties. We appreciate that the magnetic helicity evolution in the time 
interval between these epoches has to be addressed separately. In particular, 
large-scale magnetic helicity produced by galactic dynamo is antisymmet- 
ric in respect to the galactic equator while the magnetic helicity from any 
cosmological sources is obviously independent on the galactic equator po- 
sition. It is far from clear how important this asymmetry is for nonlinear 
galactic dynamos and for the observed asymmetry of magnetic field in Milky 
Way. Note also that the strong cosmological magnetic field could prevent the 
inverse MHD cascade on the scale of galaxies [15], [16]. 

We should remark that the huge helicity value (9) exists only in hot ultra- 
relativistic (T ^> m e ) early universe plasma where oi(T) is sufficiently large. 
The evolution of magnetic helicity H(t), or how cosmological magnetic he- 
licity feeds protogalactic fields is a complicated task. In the nonrelativistic 
plasma, first, positrons vanish, then with the cooling for the frozen-in mag- 
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netic field B ~ T 2 — > the electron density at the main Landau level drops, 
no- — > 0, resulting in a — > 0, and the magnetic helicity production becomes 
impossible. 

Let us note that the neutrino collision mechanism [5] can not produce 
magnetic helicity unlike in our collisionless mechanism. This immediately 
comes after the substitution of the electric field term stipulated by weak 
collisions E = — | J cxt | V/cr (a is the electric conductivity in the ultrarel- 
ativistic plasma) and taken from Eq. (5) in [5], where the electric current 
Jcxt ~ Gp is caused by the friction force due to the difference of weak cross- 
sections for neutrino scattering off electrons and positrons. This current is 
directed along the fluid velocity. The generalized momentum P = w e 7 e V, 
w e = 4T is the enthalpy, 7 e ^> 1 is the 7-factor in the ultrarelativistic plasma, 
obeys Euler equation [4] 

(d t + V • V) P = + TOt B X B + weak terms , (10) 

n e Aixn e 

from which retaining the standard MHD terms only (the first and the sec- 
ond ones in the r.h.s. of Euler equation) one can obtain the velocity V oc 
fl o [— Vp/n e + (rot B x B)/47rn e ] that does not contribute (in the lowest 
approximation over ~ G 2 F ) to the helicity change. 

Let us note that we rely here on homogeneous magnetic fields with a scale 
which is less (however comparable) than the horizon In, hence the magnetic 
force lines are closed within the integration volume f v d 3 r(...), or applying 
the Gauss theorem one can show that the contribution of the first term in 
the r.h.s. of Eq. (10) to the helicity production (6) vanishes. This is exactly 
like for the gauge transformation of the vector potential A in the helicity 
H = J v d 3 rA ■ B, A — > A + V%. On the other hand, there remains an open 
question how to define the gauge invariant helicity for superhorizon scales. 

There are other astrophysical objects for which axial vector weak forces 
acting on electric charges and driven by neutrinos can lead to the amplifica- 
tion of mean magnetic field and its helicity as given in (7). For instance, the 
neutrino flux vorticity which is proportional to Vxj(„) can vanish for isotropic 
neutrino emission from supernovas in the diffusion approximation when neu- 
trino flux j(jy)(r) is parallel to the radius r. In such case the mechanism of 
collective neutrino-plasma interactions originated by the axial vector weak 
currents becomes more efficient to amplify magnetic field than the analogous 
mechanism based on weak vector currents [17]. 
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